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The importance of taking into account inter-site f — f hybridization in electron structure calculations for 
Ce metal and cerium heavy fermion compounds was studied. We demonstrate that for heavy-fermion systems 
such as cerium compound CeCu2Si2 / — / hybridization can be neglected and Anderson model application is 
well justified. On another hand for cerium metal / — / hybridization is strong enough to provide the contri- 
bution to hybridization function comparable to hybridization between 4/ and itinerant electrons. We argue 
that in the case of Ce only the most general Hamiltonian combining Hubbard and Anderson models should 
be used. 
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The mysterious properties of metallic Ce, which has 
paramagnetic phase with the local magnetic moments at 
ambient pressure and room temperature (Ce— 7 phase), 
and show the absence of local moments and Pauli para- 
magnetism below ^100 K (Ce— a phase) rivet attention 
of the researchers [T]. For decades the electronic and 
magnetic properties of metallic Ce and heavy fermion 
cerium compounds were considered in the frameworks 
of the single impurity problem mainly using Anderson 
impurity model [2]: 



Hsiam = E £ k4<A* + £ f E PJ° + U ntfnfl + 

ker <y 

E faU* + V£fic k <r) ■ (1) 
ko- 

where localized /-electrons with on-cite Coulomb in- 
teraction term Uhffhf^, hybridize with itinerant o 
electrons described by dispersion £k with a hybridiza- 
tion strength parameter Vk- 

One can introduce noninteracting Green function Qq 
(defined as Green function with Coulomb interaction 
switched off): 



G (iuj n ) = (iu n + fi - e d - A(iuj n )) 1 



(2) 



where ui n = (2n + l)nT, n = 0, ±1,±2, ... are Mat- 
subara frequencies and hybridization function A(iuj n ) 
is defined as: 



A(iw„) = 



iu: n - e k + \i 



Then the problem that should be solved is to de- 
scribe /-electrons with on-cite Coulomb interaction in 



an effective media defined by noninteracting Green func- 
tion Qq ([2]) where interaction with effective media is de- 
termined by hybridization function A.(iu> n ) 

The calculations performed using this model allowed 
to obtain consistent description of the evolution of mag- 
netic and electronic properties as due to appearance of 
the Kondo scattering in a— phase of Ce. The impu- 
rity models were applied for the study of the magnetic 
susceptibility [3], specific heat [3] and different types 
of spectra (photoemission [3], Bremsstrahlung isochro- 
matic J5], electron-energy- loss [5]). Fitting of the theo- 
retical result obtained within impurity models to differ- 
ent experimental data (protoemission spectra, suscepti- 
bility etc.) allows to extract the most important param- 
eters in Kondo physics - Kondo temperature Tk [5]. 

While Anderson impurity model (JlJ has allowed to 
capture main energy scale in heavy-fermion physics - 
Kondo temperature Tk, it cannot describe coherence 
effects when at low temperatures rich phase diagram 
appears with long-range magnetic ordering and super- 
conductivity. Basic model used to describe such effects 
for /-systems is periodic Anderson model (PAM) with 
Hamiltonian: 

i<7 ij<7 icr 

uJ2 h n h n + E (v<Ah° + v*jfl*Zi°) ■ ( 4 ) 

It deals with localized /-electrons on all sites embedded 
in itinerant c-electrons bath with a term responsible for 
hybridization between localized and itinerant electrons. 
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Fig. Q] (Color online). LDA Ce-4/ partial DOS for Ce- 
a (dashed, red), Ce-7 (solid with circles, violet) and 
CeCu2Si2 (solid, black). The Fermi level is in zero. 

In both impurity ([T]) and periodic ^ Anderson mod- 
els hybridization between /-electrons on different lattice 
sites is assumed to be absent in contrast to Hubbard 
model where competition between inter-site / — / hy- 
bridization and Coulomb on-site interaction is explicitly 
defined: 
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If one cannot neglect inter-site / — / hybridization 
then the most general Hamiltonian combining Hubbard 
and Anderson models should be defined and studied: 

H — e c c icr Ci C + tijc\ a Cj a + 

E iiHh* + e / E 4 + u E 44 + 

ia i 

E(^4/^+^/j CT c 



(5) 



In the present paper we investigate the problem of 
applicability of Anderson model to study cerium and 
cerium compounds and estimate the strength of inter- 
site / — / hybridization. We demonstrate that while 
for heavy-fermion systems such as cerium compound 
CeCu2Si2 / — / hybridization can be neglected and An- 
derson model (U) application is well justified, for cerium 
metal inter-site f — f hybridization is strong enough 
giving contribution to hybridization function ([3]), which 
is comparable to hybridization of /-electrons with itin- 
erant electrons. In the last case only the most general 
Hamiltonian (|S|) should be used. 




Fig. [2] (Color online). Full orbital LDA band struc- 
ture for Ce-7 is shown by blue dashed curve. The band 
structure obtained by removing Ce— s,p, d states from 
the self-consistent LDA hamiltonian is shown in black. 
They are more localized, but still cannot be considered 
as atomic levels. The Fermi level corresponds to zero 
energy. 

With the use of the Linear muffin-tin orbitals 
(LMTO) method [7] and the Local density approxima- 
tion (LDA) we show that Ce-4/ states in metallic Ce 
should not be described simply as impurity levels. These 
states do form bands and /— / hopping matrix elements 
between different Ce sites are sizable. In contrast the 
/-states in Ce compounds are more localized and do not 
show significant band dispersion. 

We start from the Ce-4/ partial Density of states 
(DOS) for Ce-a, Ce-7 an d CeCu2Si2 presented in Fig.[T] 
One may see that the widths of the DOS are comparable 
for all three systems and hence it may be expected that 
the band characteristics of f-states in these compounds 
are similar. Since the similarity in the position, width 
and shape of partial DOS is most pronounced for Ce-7 
and CeCu2Si2 we will use these two systems to compare 
band effects. 

The real band structures obtained in the self- 
consistent LDA calculation for Ce-7 is shown by dashed 
curves in Fig. [5] Seven Ce-4/ are spread over wide en- 
ergy [-0.4 eV, 1 eV] (compare with Fig.[l|. 

In order to check whether Ce-4/ states can be 
treated as independent impurity states we remove (set 
zero) all the matrix elements from the self-consistent 
LDA hamiltonian except Ce-4/. The self-consistent po- 
tential for the real material is still used, so that the 
resulting band structure is not the same as for hypo- 
thetical "Ce-f-only ions" in Ce-7 type lattice. The band 
structure obtained within this method can be thought 
as the actual dispersion of Ce— / states in real Ce— 7, 
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Fig. [3] (Color online). Full orbital LDA band structure 
for CeCu2Si2 is shown by blue dashed curve. The band 
structure obtained by removing Ce—s,p,d states from 
the self-consistent LDA hamiltonian is shown in black. 
The Fermi level corresponds to zero energy. 

where hybridization with Ce— s,p, d states was switched 
off. The same procedure was previously applied for the 
analysis of the chemical bonding in Ag2Ni02 [8]. 

The comparison of full-orbital LDA bands structure 
and one obtained removing Ce-s,p, d states from the ba- 
sis set is shown in Fig. [3J One may see that the band 
dispersion of Ce-4/ states is quite similar, and that these 
states still form the real bands, rather than atomic lev- 
els. The band- width W/_ on ; y ~ 0.75 eV in Ce— 7. 

In order to show that this situation is specific to Ce 
we performed the same calculations for CeCu2Si2- The 
results are presented in Fig. [3] In contrast to the case 
of metallic Ce the absence of the hybridization between 
Ce— 4/ and Ce-s,p,d states leads to the loss of band 
dispersion. The reason for such a different behavior of 
metallic Ce and CeCu2Si2 is rather obvious: in the last 
case Ce ions are separated by Cu and Si, direct / — / 
hopping and corresponding effective bandwidth is small 
(Wf- niy ~ 0.1 eV) and the bands are dispersionless 
like atomic levels. However, the presence of sizable band 
dispersion for metallic Ce was not taken into account in 
previous model calculations. 

The value of the Ce / — / hopping parame- 
ters estimated from the band-width and tight-binding 
parametrization or more sophisticated Wannier projec- 
tion procedure [9] results in i// ~30 meV. The presence 
of small, but finite / — / hopping may lead to a number 
of consequences. The most obvious is a direct antiferro- 
magnetic exchange interaction between Ce ions propor- 
tional to 2tff/U. Together with indirect Ruderman- 
Kittel-Kasuya-Yosida (RKKY) exchange this interac- 
tion will act against formation of a coherent state. 



The most direct investigation of the effects con- 
nected with the presence of finite / — / hoppings in pure 
Ce can be performed by a numerical solution of ([5]). us- 
ing for instance Dynamical mean-held theory (DMFT) 
or its cluster extension [10]. However, already on the 
LDA level one may show that these effects should be 
important. In order to demonstrate it the hybridization 
function on the real energy axis was constructed with 
and without / — / hopping. On the first step of this pro- 
cedure one obtains Hamiltonian in the basis of Wannier 
functions in real space as described in ref. [5] . Then one 
makes zero corresponding off-diagonal matrix elements, 
perform back Fourier transform to reciprocal space, cal- 
culate density of states using tetrahedron method and 
construct hybridization function A(e) using formalism 
developed in ref. |11) : 

where p v {e) is a partial DOS, and v - orbital index. To 
avoid numerical errors partial DOS were normalized on 
unity before apply 

The plot of the hybridization function obtained in 
this way in comparison with A(e) from conventional 
LDA calculation is presented in Fig.|U One may see that 
the most significant changes in frequency are observed 
near the Fermi level. The full description of the elec- 
tronic properties of the system with given hybridization 
can be obtained only by numerical solution of many- 
body problem. However, already on the LDA level we 
obtain that the ratio Alda(£)/ At ff =o(s) is of order 2 
for Ce and 1.2 for CeCu2Si2, where A(e) is averaged 
over the region of 1 eV around the Fermi level total hy- 
bridization function. This demonstrates an importance 
of the account of direct / — / hopping matrix elements 
in a real many-body calculation 

To sum up, in the present paper we've shown that 
there is sizable f — f hopping matrix element in the 
metallic Ce. This implies that the full description of 
the electronic properties of Ce should be obtained not 
within the frameworks of the single impurity, but rather 
in lattice models, where hopping parameters between 
different f-sites are implicitly taken into account. Thus, 
multi-band Hubbard model is one of the models suitable 
for such an investigation. 
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Fig.[4](Color online). Total hybridization function A(e) 
for Ce— 7 as defined in ((SJl calculated in conventional 
LDA and in LDA, where in self-consistent Hamiltonian 
off-diagonal inter-site f — f matrix elements were put 
to zero. The Fermi level corresponds to zero energy. 
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